Quantitative chemical mapping at the atomic scale 
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Atomic-scale mapping of the chemical elements in materials is now possible using aberration- 
corrected electron microscopes but delocalization and multiple scattering can confound image in- 
terpretation. Here we report atomic-resolution measurements with the elastic and inelastic signals 
acquired on an absolute scale. By including dynamical scattering in both the elastic and inelastic 
channels we obtain quantitative agreement between theory and experiment. Our results enable a 
close scrutiny of the inelastic scattering physics and demonstrate the possibility of element-specific 
atom counting. 



PACS numbers: 

The short de Broglie wavelength of high-energy elec- 
trons enables a scanning transmission electron micro- 
scope (STEM) to form an electron beam of sub-atomic 
dimensions. Such beams can be used to probe the atomic 
and electronic structure of materials with extremely high 
selectivity. Collection of the inelastic scattering resulting 
from the excitation of atomic core electrons, via electron 
energy- loss spectroscopy (EELS), provides both chem- 
ical and electronic bonding information. Recently it 
has become possible to utilize the core-loss EELS sig- 
nal in a STEM to map the chemical composition and 
electronic structure of materials at the atomic scale [1- 
6]. This capacity constitutes an extremely powerful tool 
in materials characterization. In the present work, we 
consider whether it is possible to perform quantitative 
chemical mapping at the atomic scale. Such a capability 
would permit, for example, element- specific column-by- 
column atom counting, with significant implications for 
the materials and condensed matter communities. To 
this end, it is crucial to establish a suflficiently accurate 
interpretation of the contrast in atomic-resolution chem- 
ical maps, and, in particular, to determine whether the 
dominant inelastic scattering is treated correctly in cur- 
rent theories based on a single-particle description of the 
core- excitation process. Here we report the first atomic- 
resolution chemical mapping with the EELS signal ac- 
quired on an absolute scale (i.e. normalized with respect 
to the incident beam). We demonstrate that quantifi- 
cation within a single-particle picture is indeed possible 
under favorable circumstances. In other cases, we find 
discrepancies between the simulations and experiments 
indicative of a break-down of the single-particle model. 

In a STEM, images are generated by scanning a conver- 
gent electron beam across an electron-transparent speci- 
men in a raster- like fashion, while collecting various scat- 
tered signals in synchronization with the scanning beam. 
A common mode of atomic-resolution imaging in the 
STEM is annular dark-field (ADF) imaging, whereby the 



electrons elastically scattered to high angles are collected 
by an annular-shaped detector. While atomic-resolution 
ADF imaging of crystals can now be achieved routinely, 
quantification of the image contrast has a long and con- 
troversial history [7, 8]. Discrepancies between experi- 
mental and simulated contrast levels, referred to as the 
"Stobbs factor", were often found to be as large as 3- 
5 times, and led to speculation that fundamental scat- 
tering mechanisms were missing from the simulations. 
With improved quantification techniques, however, re- 
cent work [9] has persuasively argued that discrepancies 
in ADF image contrast can be accommodated by includ- 
ing a physically-reasonable source size in the simulations. 
A major conclusion of that work is that multislice calcu- 
lations [10] incorporating a frozen-phonon model of ther- 
mal diffuse scattering [11] correctly describe the domi- 
nant scattering mechanisms in ADF imaging. 

On the other hand, core-loss EELS oflFers a wealth of 
information not accessible by ADF imaging, including 
unique identification of the chemical elements and bond- 
ing information. However, the inelastic scattering pro- 
cesses relevant to core-loss EELS add a level of complex- 
ity over and above that of ADF imaging. In this study, we 
have performed carefully calibrated experiments, allow- 
ing atomic-resolution chemical signals to be extracted on 
an absolute scale. By simultaneously recording the elas- 
tic signal, also on an absolute scale, and using it to char- 
acterize the electron beam and sample thickness, we elim- 
inate all free parameters in the acquisition of the core- loss 
signals. Coupled with simulations that incorporate both 
core-loss inelastic scattering and dynamical elastic scat- 
tering, the present work enables a close scrutiny of the 
scattering physics in this inelastic channel. 

Experiments were performed using an aberration- 
corrected Nion UltraSTEM operating at 100 kV with a 
beam convergence semi-angle of 31.8 ±0.1 mrad. The 
specimen consisted of a DyScOs single crystal in the 
[101] crystallographic orientation. The EELS signal was 
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recorded using a Gatan Enfina spectrometer with a col- 
lection semi-angle of 79.8 ±0.1 mrad. During all ex- 
periments, ADF images were recorded simultaneously 
with the energy-loss spectra. The inner and outer ADF 
collection semi- angles were 98.2 ± 0.2 and 294.5 ±0.5 
mrad, respectively. Both the EELS and ADF signals 
were recorded so as to allow their calibration in terms of 
the fractional beam intensity. Data was obtained for a 
large range of specimen thicknesses. This permits a direct 
comparison with theoretical predictions of absolute cross 
sections and chemical maps from the Dy-M4^5 {3d 4/), 
Dy-A^4,5 (4(i Af and continuum), and Sc-L2,3 (2p 3d 
and continuum) edges. The different symmetries of the 
Dy and Sc sublattices, and the presence of two different 
Dy edges, allow us to identify and separate elastic and 
inelastic effects. 

As a first step, we performed quantification of the ADF 
images. Not only does this provide a check on various 
experimental parameters, it also enables us to determine 
the effective source size of the instrument independently 
from the EELS measurements. Since source size effects 
are manifest in the same way in both ADF and EELS 
[12], this step is crucial for the quantitative analysis of the 
core- loss signal. ADF quantification is shown in Fig. 1. 
The inset of Fig. 1(a) shows an ADF image of one unit 
cell of [101] DyScOs. The ADF image consists of bright 
spots at the positions of Dy atomic columns and weaker 
spots at the positions of Sc atomic columns (O atoms are 
not discernible). The graph in Fig. 1(a) compares the 
experimental posit ion- averaged, or mean, ADF signal as 
a function of sample thickness with the prediction from 
frozen phonon multislice simulations. The general agree- 
ment is very good (to within about 5%). A quantitative 
comparison of the contrast in experimental and simulated 
images enables us to extract the effective source size of 
the instrument. Fig. 1(b) compares the experimental and 
simulated contrasts of Dy columns, Sc columns and the 
background. Here the contrast is defined as the ADF sig- 
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FIG. 1: Quantitative determination of the effective source dis- 
tribution from the simultaneously-acquired ADF-STEM im- 
ages, (a) Comparison of the mean ADF signal with the frozen 
phonon simulation on an absolute scale, (b) Comparison of 
the experimental image contrast with the frozen-phonon sim- 
ulation incorporating the the determined source distribution. 



nal at the relevant position divided by mean ADF signal, 
eliminating any errors in the absolute cross sections. The 
effective source distribution was modeled as a combina- 
tion of Gaussian and truncated-Lorentzian distributions, 
fitted using nonlinear least-square optimization. The ef- 
fective source of our instrument is described by a Gaus- 
sian with a full- width at half-maximum (FWHM) of 0.71 
A convolved with a truncated Lorentizian with a FWHM 
of 0.18 A that reflects residual aberrations in the instru- 
ment alignment. 

The simulations used for quantitative analysis of the 
core-loss EELS signal were based on a relativistic multi- 
pole core- loss theory [13, 14], performed on a GPU clus- 
ter [15]. The relevant scattering physics is described 
by the following expression for the wave function of an 
inelastically-scattered electron at the specimen exit sur- 
face: 

M^)= J d2x"dVG„(x,x")cT„K.(x")Go(x",x')Vo(x'). 

(1) 

In this expression, ipo is the wave function of the high- 
energy electron beam at the specimen entrance surface, 
Va is a quantum-mechanical matrix element describing 
the effect of core excitation on the high-energy elec- 
tron, a a is the interaction constant, and Go and Ga 
are the Green's functions describing the dynamical elas- 
tic scattering, or "channeling", of the high-energy elec- 
tron before and after the core excitation event. Chan- 
neling was calculated using a frozen-phonon multislice 
method, while the matrix elements Va were calculated 
within a single-particle model where the initial and fi- 
nal states of the core electron are described by atomic 
wave functions [13, 16]. Since Eq. (1) includes channel- 
ing before and after core excitation, it is referred to as 
a "double-channeling" theory [17-19]. We find that al- 
though a single-channeling approximation, obtained by 
neglecting channeling after core excitation [2, 20, 21], 
produces qualitatively similar results, double-channeling 
effects are important for quantitative agreement, espe- 
cially in the presence of heavy elements such as Dy. 

Fig. 2 shows a collage of experimental and simulated 
chemical maps of [101] DyScOa for a range of specimen 
thicknesses. The majority of the maps are seen to exhibit 
local maxima at the sites of the corresponding atomic 
columns for all thicknesses in the recorded range. This is 
a highly desirable feature, as it allows an intuitive inter- 
pretation of the maps in terms of the locations of atomic 
sites. These features are well-reproduced by the simula- 
tions that incorporate the effective source distribution. 
On the other hand, simulations that do not incorporate 
the effective source (not shown) reveal that in the Dy-Ms 
and Dy-A^4^5 maps local minima do, in fact, develop at 
the Dy columns as the thickness increases beyond about 
20 nm. While similar minima have been observed previ- 
ously and were attributable to coherent effects [2, 19, 22- 
26], our large collection angle implies that here the min- 
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FIG. 2: Quantitative atomic-resolution chemical mapping of [101] DyScOs. For each chemical signal (indicated atop) and 
specimen thickness (indicated left) we show raw, noise- filtered, unit- cell- averaged, and simulated maps. Unit-cell-averaged 
maps are from about 25 unit cells. Intensity scales refer to the fractional beam intensity. The experimental spectra and 
energy-loss integration windows are shown in Fig. 3(a). Simulations incorporate the overall scaling factors from Fig. 3(b). 



ima develop as a consequence of strong elastic scattering 
by the heavy Dy column prior to core excitation, with 
the scattering subsequent to core excitation also playing 
a role. For the Dy-A^4 5 maps, the form of the matrix ele- 
ments is such that the minima are even more pronounced, 
to the point where the Dy-7V4 5 maps in Fig. 2 can ex- 
hibit minima at the Dy columns even when the effective 
source is included. The effect is particularly strong in the 
Dy-A^4 5 on-resonance map, which exhibits clear minima 
across almost the entire range of thicknesses. Contrary 
to this behavior, the Sc-L2,3 maps exhibit no such min- 
ima because the lighter Sc columns imply weaker elastic 
scattering. Overall, the results in Fig. 2 demonstrate very 
good agreement between the simulated and experimental 
chemical maps across the range of recorded thicknesses. 

Let us proceed to a quantitative comparison of the ex- 
perimental and simulated chemical maps. The energy- 
loss spectra and integration windows of the signals are 
shown in Fig. 3(a). For Dy-Ms and Sc-L2,3 we have 
analyzed the near-edge signals since they are strongest, 
while for Dy-A^4^5 we have avoided complications asso- 
ciated with the Fano resonance and have analyzed the 
signal at 60 eV above threshold. We proceed by com- 
paring the absolute cross-sections, as represented by the 
mean inelastic signals as a function of sample thickness, 
as shown in Fig. 3(b). Firstly, it is noted that for Dy- 
M5, the simulation significantly overestimates the mean 
signal (dashed line in Fig. 3(b)). The discrepancy is at- 



tributable to inelastic scattering associated with multi- 
ple plasmon excitations, which is absent from the sim- 
ulations. Ideally, the effects of such scattering would 
be removed from the experimental spectra by deconvolu- 
tion. However, we find that this method is unreliable in 
the case of atomic-resolution EELS where the low signal 
means that the spectra are noisy and the extended edges 
are not fully recorded. Instead, we chose to incorporate 
the effects of multiple plasmon scattering into the simu- 
lations as a correction given by Beer's law e~^l^^ ^ where 
i is the sample thickness and Aq is the inelastic mean-free 
path. This correction has been applied to the simulated 
Dy-Ms and Sc-L2,3 signals in Fig. 3(b). After applying 
the correction, we observe that the thickness-dependent 
trends of the simulated mean signals agree well with the 
experiments. Moreover, the absolute value of the Dy-Ms 
mean signal is in agreement with experiment to within 
about 5%. For the Sc-L2,3 signal, we observe that the 
simulation overestimates the signal by about 37%. For 
the Dy-A^4^5 off-resonance signal, no Beer's law correc- 
tion was applied because the integration window receives 
signal from the resonance while losing signal at the same 
time. The simulation overestimates the Dy-7V4^5 signal 
by a factor slightly less than 3. As discussed below, the 
larger discrepancies for Sc-L2,3 and Dy-A^4^5 indicate a 
break-down of either the single-particle model or assump- 
tions in the processing of the experimental data. Overall, 
the mean inelastic signals are in agreement with the sim- 
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FIG. 3: Quantitative comparison of the inelastic signal and 
chemical mapping contrast, (a) Background-subtracted ex- 
perimental spectra at the Dy-Ms, Sc-L2,3 and Dy-A^4,5 edges, 
(b) Comparison of experimental and simulated mean inelas- 
tic signals, (c) Comparison of experimental and simulated 
chemical mapping contrast. 



ulations within the expected uncertainties, but highlight 
the need to include multiple inelastic scattering. 

Similar trends of agreement carry over to the contrasts 
in the chemical maps, as presented in Fig. 3(c). For ex- 
ample, the Dy-Ms contrast is reproduced closely by the 
simulation (when the latter includes the effective source 
distribution). The Sc-L2,3 contrast is slightly overesti- 
mated by the simulation, the discrepancy being around 
25% for thicknesses below 50 nm, and reducing to about 
15% for thicknesses greater than 50 nm. For the Dy- 
A^4^5 off-resonance maps, the discrepancies are more sig- 
nificant, with the experimental contrast being lower by 
nearly a factor of 2. Nonetheless, the overall thickness- 
dependent trends of the experimental and simulated con- 
trasts are in agreement for both Sc-L2,3 and Dy-A^4^5. 

The discrepancies observed for Sc-L2,3 are likely due 
to one of three factors. Firstly, the ^2,3 branching-ratio 
anomaly of Sc indicates that there is a large core hole- 
3d interaction, implying that the single-particle model 
used in the simulations is not fully correct. Secondly, the 
final states for the Sc-L2,3 signal are sensitive to solid 
state effects not included in the simulations. In line with 
these observations, comparison of the experimental and 
simulated Sc-L2,3 mean signal at 60 eV above thresh- 
old (not shown) exhibits substantially better agreement. 
Thirdly, the Sc-L2,3 edge is sitting on a background sig- 
nal whose spatial distribution changes with energy loss. 
In the processing of experimental data, an extrapolation 
of the pre-edge intensities is used to extract the core- 
loss signals. The assumption of this procedure is that 
the extrapolated background has the same spatial distri- 



bution as the pre-edge signal, and this assumption can 
be invalid at low energy losses. In the case of Dy-A^4^5, 
the observed discrepancies are likely caused by the strong 
Fano resonance [27-29]. As shown in Fig. 3(a), the Dy- 
^4,5 pre-edge intensities are negative due to the Fano 
resonance, making the background subtraction problem- 
atic. In addition, as exemplified in Fig. 2, the signal in 
the Dy-A^4^5 on-resonance maps spreads throughout the 
unit cell due to the relatively large spatial extent of the 
4:d and 4/ orbitals of Dy. Even though the Dy-A/4 5 off- 
resonance maps were extracted from about 50 eV above 
the resonance, it is possible that a significant fraction of 
the signal originates from the resonance. In contrast to 
the Sc-L2,3 and Dy-A^4^5 signals, the situation for Dy-Ms 
is relatively simple. For example, the atomic-like final 
states for the Dy-Ms signal fully justify the single-particle 
atomic model used to compute the quantum-mechanical 
matrix elements. Moreover, the relatively large energy 
loss gives rise to a slowly- varying background signal, en- 
suring validity of the background subtraction. 

In conclusion, we have demonstrated the simultane- 
ous acquisition of quantitative ADF images and EELS 
maps in an aberration-corrected STEM. Coupled with 
double-channeling simulations of core-loss inelastic scat- 
tering, this has enabled a quantitative examination of 
the core-loss signal utilized in atomic-resolution chemical 
mapping. We found that by taking into account the ef- 
fective source distribution determined from the ADF im- 
ages, both the absolute signal and the contrast in atomic- 
resolution Dy-Ms maps can be closely reproduced by the 
double-channeling simulations. At lower energy losses, 
discrepancies are present in the Sc- 1/2,3 ^ind Dy-A/4^5 
maps due to the energy-dependent spatial distribution of 
the background spectrum, core-hole effects, and omitted 
complexities in the final states. This work has demon- 
strated the possibility of using quantitative STEM-EELS 
for element-specific column-by-column atom counting at 
higher energy losses and for atomic-like final states, and 
has elucidated several possible improvements for future 
work. 
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